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Abstract 

We construct a gauge-fixing procedure in the path integral for gravitational 
models with branes and boundaries. This procedure incorporates a set of gauge 
conditions which gauge away effectively decoupled diffeomorphisms acting in the 
(d + l)-dimensional bulk and on the <i-dimensional brane. The corresponding 
gauge-fixing factor in the path integral factorizes as a product of the bulk and 
brane (surface-theory) factors. This factorization underlies a special bulk wave- 
function representation of the brane effective action. We develop the semiclas- 
sical expansion for this action and explicitly derive it in the one-loop approxi- 
mation. The one-loop brane effective action can be decomposed into the sum of 
the gauge-fixed bulk contribution and the contribution of the pseudodifferential 
operator of the brane-to-brane propagation of quantum gravitational perturba- 
tions. The gauge dependence of these contributions is analyzed by the method 
of Ward identities. By the recently suggested method of the Neumann-Dirichlet 
reduction the bulk propagator in the semiclassical expansion is converted to the 
Dirichlet boundary conditions preferable from the calculational viewpoint. 



1. Introduction 

In this paper we extend the method of quantum effective action in brane models £Q 
to gravitational systems invariant under general coordinate diffeomorphisms. Specifi- 
cally, we will be interested in the peculiarities of the gauge-fixing procedure caused by 
the presence of branes/boundaries. It is well known that branes break the full diffeo- 
morphism invariance and give rise to the dynamical brane bending modes [2] which 
produce ghost instabilities [3J E] and generate a low strong-coupling scale 0. This 
leads to difficulties in the construction of long-distance modifications of gravity theory 
that could underly the dark energy phenomenon both at the classical and quantum 
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level |H] . Quantum effects in gravitational and, in particular, cosmological brane mod- 
els might also be important for the mechanism of the cosmological acceleration that 
can be facilitated by the 4D conformal anomaly in the presence of extra dimensions [7j. 
So here we want to study the gauge-fixing aspects of such problems in the quantum 
domain. 

Peculiarities of the gauge- fixing procedure in brane models follow from the structure 
of their action 

S[G AB (X)} = S B [G AB (X)} + S h [g al3 (x)}. (1.1) 

It contains the bulk and brane parts as functionals of the bulk metric Gab{X) and 
the induced metric on the brane g a p{x). The bulk (d + l)-dimensional and the brane 
d- dimensional coordinates are labeled respectively by X = X A , A = 0, 1, ...d, and 
x = x a , a = 0,1,. ..d — 1, the brane is embedded into the bulk by means of some 
embedding function e(x), 

X A = e A (x), (1.2) 

and the induced metric reads in terms of the bulk metric as 

9c#{x) = d a e A (x) G AB (e(x)) d p e B {x). (1.3) 

Generically, the bulk part in (|1.1|) is the Einstein action with the cosmological term 

S B [ G AB (X) ] = | jT d d+l X G 1 ' 2 (R(G) - 2A) - 2 jf d d xg x l 2 (1.4) 

which includes the integral over the bulk B and the Gibbons-Hawking integral over the 
brane b. The latter involves the trace of the brane extrinsic curvature K = g a,3 K a p 
and is necessary for the consistency of the variational procedure for this action. 1 

The brane part of (jl.4|) is given by a covariant ci-dimensional integral over b. De- 
pending on the model it contains the brane tension term (in the Randall-Sundrum 
model |Sj), the brane Einstein term (like in the Dvali-Gabadadze-Porrati model [Hj), 
the combination thereof ^U] or other covariant structures in brane and bulk curva- 
tures Both bulk and brane parts of the full action may also contain matter fields 
which for brevity we do not consider here. Without loss of generality they can be 
included into the sets of bulk and brane metric coefficients. Also, it is worth noting 
that Eqs. (jl.l|) - (|1.4|) equally well describe the system in the bulk domain B with the 
boundary b = <9B and the brane system with the bulk satisfying the ^-symmetry with 

The extrinsic curvature is defined as a projection onto the brane of the tensor Kab = V ' atib with 
the inward pointing vector ub normal to the brane. 
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respect to the brane b. This formalism also applies in context of Euclidean quantum 
gravity with a closed compact boundary b. 

The bulk and brane parts of the action (jl.lj) are invariant with respect to diffeomor- 
phisms generated by the bulk E A (X) and brane vector fields. The corresponding 

transformations have the form 



where denotes the covariant derivative in the bulk and D a is a covariant derivative 
on the brane. The bulk diffeomorphism ()1.5j) preserves the bulk action only when 
the projection of normal to the brane, is vanishing. In contrast, the 

tangential projection of H Q (X), on the brane can be arbitrary. Because the 

bulk and induced brane metrics are not independent, this projection generates the 
brane diffeomorphism and coincides with the vector field £, a (x). Thus, the boundary 
conditions for the above gauge transformations read 



Here and in what follows the vertical bar indicates that the function in the bulk, 
4>{X), is restricted to the brane and labeled by the corresponding low case letter x, 
4>{X) \ = <f)(e(x)), with the aid of the embedding function (jl.2|) . 

The construction of the quantum effective action for any brane model should in- 
corporate the gauge-fixing procedure for the diffeomorphism symmetry (jl.5j) - (jl.7|) . 
Schematically, this corresponds to introducing the Feynman-DeWitt-Faddeev-Popov 
gauge-fixing factor in the path integral ^2] 



Here integration runs over the metric coefficients in the bulk and on the brane because 
the metric on the brane is not fixed and is subject to quantum fluctuations. The 
gauge-fixing factor should gauge away local diffeomorphisms by imposing a certain set 
of gauge conditions and introducing necessary ghost factors maintaining unitarity and 
gauge independence of the effective action. 

Curiously, despite a well-known form of the gauge algebra ()1.5|) - ()1.7J1 . no exhaus- 
tive formulation of gauge-fixing procedure for (|1.8|) is known in current literature. In 
particular, the ghost factors require the knowledge of boundary conditions for the cor- 
responding Faddeev-Popov operators acting on the space of gauge parameters E A (X). 
In view of ()1-7|) their tangential components are arbitrary and should be inte- 

grated over, and it would seem that this integration gives rise to the Neumann boundary 



S a G AB (X) = V A 5 B (X) + VbZa(X), 



(1.5) 
(1.6) 



Z<*(X) =r(x), E X (X) =0. 



(1.7) 




(1.8) 
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conditions for E a (X)\. However, this naive prescription turns out to be incorrect. As 
we show below, the derivation of the boundary conditions for ghosts is less straight- 
forward and much deeper relies on the details of gauging away the transformations 

In its turn, the preferable choice of gauge conditions in brane models is constrained 
by additional requirements. One requirement follows from the fact that within the 
brane concept the fields in the bulk are usually integrated out, so that the physics of 
the system is effectively probed only by the variables living on the brane. 2 Thus, a 
very important aspect of the effective action in brane theory is that it should be a 
functional of the induced metric g a p(%), rather then the full set of metric coefficients 
Gab(X). Therefore, the preferable gauge conditions should be chosen in such a way 
that the calculational technique for r[g a p(x)] is manifestly covariant with respect to 
brane diffeomorphisms (jl.6|) separately from the bulk ones (J 1.51) . 

The effective decoupling of brane diffeomorphisms from the bulk diffeomorphisms 
can be attained by imposing two sets of gauge conditions — brane d-dimensional gauges 
X^{ x ) — and the bulk (d+ 1) -dimensional gauges H A (X) = which gauge away their 
corresponding diffeomorphisms. Moreover, they can preserve the manifest covariance 
of the calculational scheme if they are chosen in the DeWitt background-covariant 
form ^2]- Under the splitting of the brane metric into the background (mean field) 
and perturbation parts, g a p — > g a p + h a p, this brane gauge reads as 

X^D'h^-^h, h = g^h, u , (1.9) 

where the covariant derivatives are defined with respect to the background metric 
g^y which also raises the indices of and D^. The corresponding Faddeev-Popov 
operator which we denote by J„(D) acts on the vector space of £ u (x). In the lowest 
order in it reads as a covariant <i-dimensional d'Alembertian modified by the Ricci- 
curvature potential term 

J»{D) = n^SH + R£, D (d) = gVDpDv. (1.10) 

A similar DeWitt bulk gauge condition arises under the splitting of the full bulk 
metric, Gab Gab + Hab, and reads 

H A = V B H BC - l -V A H, H = G AB H AB , (1.11) 

2 This means that in the classical action of the path integral (|1.8|) one must add sources conjugated 
only to the variables located at the brane, that is to the components of the induced metric g a f3(x). 
The Legendre transform with respect to these sources then leads to the effective action as a functional 
of the mean induced metric, r[ {g a p{x)) ]. As is well known, in the one-loop approximation this field 
is called the background, and the whole procedure reduces to the subtraction from S[ Gab(X)} in 
<|1.8f) the term linear in the deviation of the integration variable from this background - the stationary 
point of the path integral. 
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where Va is also defined with respect to the background metric Gab- The correspond- 
ing Faddeev-Popov operator, denoted by Q^(V) and acting on the space of E B (X), is 
given by 



This operator is of the second order in derivatives and, as we show below, it should 
be supplied by the Dirichlet boundary conditions for all components of This 
guarantees the decoupling of the bulk and brane diffeomorphisms even despite their 
entanglement via the relation (J1.7)) . 

Another requirement to the background field formalism for brane theories is the re- 
duction of their Feynman rules to those of the Green's functions with Dirichlet bound- 
ary conditions. As discussed in £Q such boundary conditions are calculationally much 
simpler than the original brane matching conditions of the generalized Neumann type. 
This reduction is possible in view of the duality relations between the Dirichlet and 
Neumann boundary value problems suggested in [15] . and as we show below this re- 
duction also works for gauge-fixed framework. In particular, the decomposition of the 
quantum effective action in the sum of the bulk and brane-to-brane contributions ob- 
tained in pQ, which incorporates this Neumann-Dirichlet reduction, works both in the 
gauge- field and ghost sectors of (jl.8|) . 

Below we develop the gauge-fixing procedure of the above type. We derive the 
gauge- fixing factor in Eq. (jl.8|) . show its bulk-brane factorization property underlying 
a special bulk wavefunction representation of the brane effective action ^3] and build 
the semiclassical expansion for this action with its Neumann-Dirichlet reduction. 

The paper is organized as follows. In Sect. 2 we list the main results of this paper. 
Sect. 3 presents the spacetime foliation associated with its brane/boundary and the cor- 
responding canonical formalism which facilitates the derivation of relevant boundary 
value problems. In Sects. 4 and 5 we derive the gauge-fixing factor of (|1.8|) and the 
associated bulk wavefunction representation of the brane effective action. In Sect. 6 we 
derive the Feynman diagrammatic technique for this action and explicitly present it in 
the one-loop approximation as a sum of special bulk and brane contributions. They are 
based on a mixed Dirichlet- Neumann boundary value problem for the graviton propa- 
gator. Sect. 7 presents the method of Ward identities which allow one to demonstrate 
for both of these contributions their manifest gauge independence. In Sect. 8 the di- 
agrammatic technique is completely reduced to the propagator with strictly Dirichlet 
boundary conditions. The concluding section summarizes the obtained results and dis- 
cusses their possible applications. Three appendices describe the properties of the local 
measure in the path integral for brane models, the Gaussian integration over the func- 
tional space with mixed Dirichlet-Neumann boundary conditions and the variational 
problem for the Green's function subject to such boundary conditions. 
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2. Main results 



The main results of this paper include the selection of bulk and brane gauge condi- 
tions, their corresponding gauge- fixing factor in the path integral ([1.8)1 . the one-loop 
expression for the brane effective action and its Neumann-Dirichlet reduction. In this 
section we briefly list these results that will be derived in the following sections. 

2.1. Bulk and brane gauge conditions 

In view of ()1-7|) a generic diffeomorphism in the bulk can be represented as a composi- 
tion of two transformations. One transformation is uniquely fixed by a given vector field 
on the brane by continuing its <i-dimensional diffeomorphism from the brane to 
the spacetime bulk. Another transformation is a generic (d+ l)-dimensional diffeomor- 
phism which reduces to the identity at the boundary, E A (X) \ = 0. Correspondingly, 
the gauge-fixing procedure can be split into two stages. At the first stage special bulk 
gauge conditions gauge away this (d + l)-dimensional diffeomorphisms vanishing at 
the brane. At the second stage brane gauge conditions completely fix the residual 
d- dimensional diffeomorphisms at the boundary. This strategy allows one to decouple 
the boundary diffeomorphisms from the bulk ones and render the whole formalism 
manifestly covariant with respect to the both types of diffeomorphisms. 

For this purpose the bulk gauge conditions should be special in the sense that they 
should not overconstrain the relevant (d + l)-dimensional transformations and leave 
the boundary values S M (A) | = arbitrary. This means that £ M (x) play the role 

of boundary conditions for zero modes of the relevant Faddeev-Popov operator — the 
transformations which leave the bulk gauge conditions invariant. If we denote these 
(d + l)-dimensional gauge conditions by 



then their gauge transformation serves as the definition of the bulk Faddeev-Popov 
operator Qb(V) 



To fix the bulk diffeomorphisms this operator should be invertible under the Dirich- 
let boundary conditions S A (A) | = 0, but it must have nontrivial zero modes Sj^(X), 
<5s(V) Eq(X) = 0, subject to inhomogeneous boundary conditions (|1.7|) "enumerated" 
by all possible brane vector fields ^(x). This is possible when the Faddeev-Popov oper- 
ator in ()2.2|) is of the second order in derivatives transversal to the spacetime boundary. 



H A (X) = H A (GcD(X) J d B G CD (X)) 



(2.1) 



6 S H A (X) = Q A (V)E B (X). 



(2.2) 
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The gauge conditions which generate the Faddeev-Popov operators with such proper- 
ties can be background-covariant from the bulk (d + l)-dimensional viewpoint and, in 
particular, given by the DeWitt gauge in the bulk 

The residual gauge transformations with the parameters £ M are gauged away by 
imposing the brane gauge conditions x^i x ) = 0. In order to decouple the boundary 
diffeomorphisms from the bulk ones, one must impose these gauge conditions only upon 
the brane metric 

X^x) = x^9a ( s(x),d ll g al3 (x)). (2.3) 

Similarly to ()2.2|) they generate the brane Faddeev-Popov operator J%(D), 

5^(x) = J^D)C(x), (2.4) 

which is now determined entirely in terms of quantities induced on the brane. A 
particular example of such gauge conditions which are background-covariant from the d- 
dimensional viewpoint is given by the DeWitt gauge conditions fjl.9|) with the Faddeev- 
Popov operator (|1.10|) which is nondegenerate under appropriate boundary conditions 
at the infinity of the brane (or in view of the closed compact nature of the boundary 
in Euclidean context when fjw becomes a Laplacian). 

2.2. The gauge-fixing factor 

For the gauge- fixing procedure (j2.1|) - (j2.4j) the gauge- fixing factor in ()1.8|) factorizes into 
the product of the corresponding bulk and brane factors 



( gauge-fixing )=S[H] Det D Q x 5(x) det J. (2.5) 

Here and S(x) denote respectively the (d + l)-dimensional and <i-dimensional 

functional delta-functions, 



xgb a 

and Detj)Q and det J are the corresponding ghost functional determinants of the bulk 
and brane Faddeev-Popov operators defined by (|2.2|) and (|2.4|) . To distinguish between 
the functional dimensionalities of these determinants we denote the determinant of the 
(d+ l)-dimensional theory by Det = Det^ d+1 ' ) and that of the c?-dimensional theory by 
det = Det (d) . 
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The subscript in DetnQ indicates that the functional determinant of the second- 
order differential operator Qg(V) is taken subject to the Dirichlet boundary conditions 
on the brane. The variational definition of such a determinant in terms of the Dirichlet 
Green's function of Qb(V) is given below, see Eq. ()4.9|) . In contrast, boundary condi- 
tions for J(D) in det J are not important for us, because the brane boundary is either 
absent, as in Euclidean context with a closed compact boundary of a bulk domain 
(the boundary of the boundary), or lies at infinity where the boundary conditions are 
trivial. 

One can use the t'Hooft's method of transition from the degenerate (delta-function 
type) gauges to the gauge-breaking terms both in the bulk sector 

6[H] ^(Detc AB ) 1/2 exp(^ g B b [G]), (2.8) 

SliG(X)] = -1 JjXH A (X)c AB H B (X), (2.9) 

and brane sector 

6( X ) (detc^) 1/2 exp(z5 g b b [(7]), (2.10) 
S* h [9(x)] = -\ J dxx^c^x^x). (2.11) 

Here cab and c^ u are respectively the bulk and brane gauge-fixing matrices. We use 
for them the same notation differing only by the type of indices (AB vs fiu), which 
should not lead to a confusion in what follows. 

The factorization of the gauge-fixing factor ()2.5|) underlies a special bulk wavefunc- 
tion representation of the brane effective action, see Eqs. (|5.2j) - (|5.3)l below. 

2.3. One-loop brane effective action 

The semiclassical expansion for the brane effective action within the gauge-fixing pro- 
cedure of the above type has the form 

r[ 9aP {x) } = S B [G%[ g(x) } ] + S h [g a p(x) } 

+hr 1 _ loop [g aP (x)] + 0(h 2 ). (2.12) 

Its tree-level part (the first line of this equation) follows from the classical action (|l.lj) 
calculated on the solution of classical equations of motion, G AB [g(x)}. This solution 
satisfies the bulk gauge (j2.1jl and is subject to the boundary condition on the brane - 
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the induced <i-dimensional metric g a g(x), 

6SB W = o, (2.13) 
8Gab{X) g=g o 

H A (G°(X),dG°(X)) = 0, (2.14) 

G%{X)\ h = g aP {x). (2.15) 

The one-loop part of the action is built in terms of the inverse propagator of the 
theory in the bulk, which is non-degenerate due to the contribution of the bulk gauge- 
breaking term 

f .^ m(V)j (X,X^ j;g» xr ,2,6) 

The Green's function of this operator, which determines the bulk propagator, 

Gdn(X,X') = G D ^ BiCD (X,X f ), (2.17) 

satisfies the following problem with the set of mixed boundary conditions of the gen- 
eralized Dirichlet-Neumann type 

F(V) G DN (X, X') = 5(X, X'), (2.18) 
G D a l CD (X,X')\ x = J (2.19) 

H(V)G DN (X,X')\ x = 0. (2.20) 

Here (I2.19J1 means that the induced metric components of the Green's function (a/3- 
components among A-B-components) vanish on the brane, whereas the rest of the 
boundary conditions imply the vanishing of the linearized gauge conditions. The dif- 
ferential operator of linearized gauge conditions, acting upon the Green's function in 
(12~2UD . H E ' AB (V)G I f B)CD (X,X'), is defined by 

H(v)s(x,x f ) = h e > ab (V)5(x,x>) = ™2{x>y (2 ' 21) 

In terms of these quantities the one-loop effective action reads as the sum of the 
bulk and brane effective actions both including their relevant ghost contributions (cor- 
responding to the factorization of the gauge-fixing factor (|2.5|) ) 



iA_i 00 p [ 9ap{x) ] = - r Tr DN In F + Tr D In Q 



G=G (g) 



— — tr In F DN + tr In J. (2.22) 
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Here Tro and Ttdn denote the functional traces of the bulk theory subject to the 
Dirichlet and mixed Dirichlet-Neumann (cf. (|2.18j) - (|2.2Uj) ) boundary conditions, while 
tr is a functional trace in the boundary <i-dimensional theory. _F DN is the brane-to- 
brane operator 



F w = F^ S (x,x') 

( WG DN W\\) aP,lS (x,x') + K a ^ s (D)5(x,x'). (2.23) 



This is a gauge theory generalization of the brane operator introduced in pQ. The 
first term here implies that the kernel of the Dirichlet-Neumann Green's function is 
being acted by the operators W(V) upon both arguments in the directions indicated by 
arrows. The double vertical bar indicates that both points of the kernel are restricted to 
the brane and labeled by the corresponding low case letters. That is, if the embedding 
of the boundary /brane in the bulk is denoted by X = e(x), then this explicitly means: 

WG DN W\\ 0,x') = W(y)Gvx(X,X')W(V) . (2.24) 

X=e(x), X'=e(x') 

The first order differential operator W(V) is the Wronskian operator 

W(V) = W AB ' CD (V), (2.25) 

which is determined for the symmetric second-order differential operator (J2.16|) by the 
following Wronskian relation valid for two arbitrary test functions ipi^{X) 

J dX (0i F{ V)0 2 - 0i F{ V) h)=~ J dx (0i W{ V)0 2 - 0i W(V) 2 ) • (2.26) 

B dB 

With these definitions the first term of (J2.23)) is given by the (a/3, 7<5)-block of 

the matrix (W G D n W \\) ab ' cd . The second term of (l2~2H is a contribution of the 
brane ci-dimensional part of the classical action (jl.lj) and the brane gauge-breaking 
term (gTTJ 

*#>*{p) 5(x, x') = ^ 2 ^ b + ^) (2 . 27) 
The gauge-breaking term here makes the whole operator (|2.23j) nondegenerate. 



2.4. Neumann-Dirichlet reduction 

The equation (J2.22)) generalizes the technique of to brane gravitational models 
invariant with respect to local diffeomorphisms. The main goal of £Q was a complete 
reduction of the calculational technique from Neumann-type boundary conditions to 



10 



much simpler Dirichlet ones. However, in the algorithm ()2.22j) this goal is not yet 
achieved, because of the complexity of boundary conditions for G DN (X, X') in (|2.18j) - 
(J2.20|) . Such a reduction can be done, and in the one- loop approximation it is given by 
the expression alternative to ()2.22|) 

«A-ioo P ] = —z Tr D InF + Tr D InQ 

--tr lnF D + tr In J. (2.28) 

Here the functional traces of all bulk operators are calculated subject to Dirichlet 
conditions, and the brane-to-brane operator F D is given by the expression similar to 
()2.23j) but with the Dirichlet Green's function Gd replacing Gdn 

Fjf 75 (x, x') = - ( W G D W 1 1 ) aP ' yS (x, x') + K af3 > ^ (D) 5(x, x') . (2.29) 



3. Canonical formalism in the brane foliation of 
spacetime 

The derivation of the above results is much easier in terms of the DeWitt condensed 
notations [T2|. In these notations the bulk metric is labeled by the condensed index 
a = (AB, X) including both tensor labels and the bulk coordinates X 

G a = G AB (X), a = (AB,X). (3.1) 

If we apply the same convention to the definition of the vector field S A = E A (X), 
A = (A, X), then the gauge invariance of the bulk action can be written down in terms 
of the generators of diffeomorphisms R a A 

5 s G a = R a A E A , (3.2) 

R/ = 0, (3.3) 
5G a v ' 

where the contraction of condensed indices implies also integration over the bulk coor- 
dinates X. In order to display the differential structure of these generators for various 
metric components 

R a c = R A B,c(V)5(X,X'), R AB , C (V) = 2G {AC V B) , (3.4) 

we will need the canonical formalism of the bulk gravitational action, associated with 
the brane-type foliation of the full spacetime. In this formalism the role of time is 
played by the extra-dimensional coordinate y. In view of the spacelike nature of y this 
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variable has nothing to do with a real dynamical evolution, but the local differential 
properties of the diffeomorphism transformations are very similar to those of the usual 
canonical formalism in the physical time t = X°. 

So we assume that the bulk is foliated by the surfaces of constant y, whose em- 
bedding is determined by the embedding functions, X A = e A (x a ,y), including the 
embedding of the brane (|1.2j) . e A (x a ) = e (x a ,0). Such a foliation determines the 
vielbein of vectors e A and n A respectively tangential and normal to the slices, e A = 
d a e A , nA_e A = 0, GABn A n B = 1. It also gives rise to the induced metric on the 
brane and other surfaces of constant y - the generalization of Eq. (jl.2|) . g a p{x,y) = 
y) G AB (e(x, y)) ef(x, y). 

This foliation generates the lapse and shift functions defined as normal and tangen- 
tial projections of the local "velocity" vector d y e (x, y) = e A with which the bulk slice 
evolves in "time" y 

N ± = n A d y e A (x,y), N a = N p , N p = ef G BA d y e A {x, y). (3.5) 

In fact these functions are equivalent to G^-components of the full metric, so that 
Gab(X) can be parameterized in terms of g a p and N A = (iV", iV -1 ). Then the bulk 
Einstein action (jl.4|) has a well-known ADM form in terms of the extrinsic curvature 
K a p of constant y slices and their scalar curvature (g) 

SB[G]= l^G J dy J dx 9 1/2 N ± {K 2 -K^ + R^(g)-2A), (3.6) 

y>o 

K a/3 = ( dygap - D a Nf3 - D p N a ) . (3.7) 
It is important that only the "velocities" of g a p, 

dyda/3 = 9a/3, (3.8) 

enter the Lagrangian, while the lapse and shift functions N A are not dynamical and 
serve as Lagrange multipliers in the y-time canonical formalism of the action ()3.6j) . In 
spacetime condensed notations of (J3.1j) we will denote the induced metric of y-slices by 
g l and this decomposition will look like 

G a = (g\N A ), a=(i,A), i={a(3,X), A=(a,±,X), 

g* = gaP (X), N A = (N a (X), N^X)). (3.9) 

In what follows we will also need canonical condensed notations in which the indices 
(a, i, A) include together with discrete labels only the brane coordinates x, and the 
contraction of these indices implies the integration only over x. 

G a (y) = (g\y),N A (y)), i = (a/3,x), A = (a, _L, x), 

g l (y)=g a p(x,y), N A (y) = (N a (x,y), N ± (x,y)) (3.10) 
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In these notations the generators of Eq. (j3.2)l . R A , form delta-function type kernels in 
the variable y with two entries a — > (a, y), A — > (A, ?/), 

i^ = i2 B x(^)*(y-lO- (3- 11 ) 
Various components of R a A (d y ) are either the ultralocal (multiplication) or differential 
operators acting on the first argument of the delta function. 

For sake of brevity, when using the condensed notations of the canonical or space- 
time nature we will not introduce special labels to distinguish between them. As a rule, 
when the y-argument is explicitly written we imply that the corresponding condensed 
indices are canonical, i.e. they contain only discrete labels and brane coordinates x, and 
their contraction does not involve implicit y-integration. For example, the left-hand 
side of ()3.3j) can be written down in the form 

«/|g=™^ A^(A,y), (3.12) 

where the integration over y (implicit in the contraction of the spacetime condensed 
index a) removed the delta function contained in R A a and the result boiled down to 
the action of the differential operator R A a (d y ) on 5S/SG a (y). This operator obviously 
differs from that of Eqs. ()3.2|) and (J3.11|) by the functional transpositon - integration 
by parts, because in contrast to ()3.2|) it acts on the test function with respect to the 
upper index a. This fact is indicated by the order of operator indices reversed relative 

to Eq. (mm). 

Another distinction between these two types of condensed notations concerns func- 
tional derivatives. We shall always reserve the functional variational notation 5/5G a = 
5/5G a (y) for the variational derivative with respect to the functions of y, while the vari- 
ational derivative with respect to the functions of brane coordinates will be denoted 
by partial derivatives. For example, S/SG a = 5/5G A b(X) vs d/dg l = 5/5g a p{x). 

In these notations the action has the form 

S B [G] = J dyL B (g,g,N) (3.13) 

where the ADM Lagrangian is displayed explicitly depending on g % — G l , and all x- 
derivatives are implicit in condensed canonical notations. With the definition of the 
momenta conjugated to g % 

Pi = Pi{g,g,N) = — (3.14) 

the bulk action can be rewritten in the canonical form 

S B [g,p,N} = J dy{p t g l -N A H A (g,p)), (3.15) 
S*[g, P ,N]\ p=p0{g . >N) = S*[Gl (3.16) 
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where HA{g,p) = (H a (g,p), H±(g,p)) is a set of momentum and Hamiltonian con- 
straints. Similarly to the ADM formalism in the physical time, these constraints as 
functions on the phase space of {g l ,Pi) satisfy the Poisson bracket relations 



{H A ,H B } = U^ B H C (3.17) 

with structure functions U AB = U^ B (g), signifying that these constraints belong to the 
first class according to the Dirac classification. As a consequence of (|3.16|) they also 
comprise the (_L ^-projections of the Einstein equations in the bulk, 

^ , \ i 5S B \g,N] 

H A (g,p) \ p=p0{g>m = jj^t- ■ (3.18) 

Due to the constraint algebra the canonical action (J3.15j) is invariant under the 
gauge transformations with local (arbitrary time and space dependent) parameters 
E A (X) satisfying S ± (X)| = 0. 3 These transformations are canonical and, therefore, 
ultralocal in y for phase space variables, but involve the y-derivative of S A (X) for 
Lagrange multipliers ^H] 



5 E g* = {g\H A }~ A , 5 E p t = { Pl , H A } ~ A , (3.19) 

6 B N A = ^A_ jj^ N B S C ^ 20 ) 

In view of the relation (j3.16j) between the canonical and Lagrangian formalisms 
various components of R a A follow from the transformations ()3.19j) - (j3.20j) [TK] 



r) / / 



(3.21) 

P=P°(g,g,N) 



dp 

R B a = R B A (dy) KV - y') = (8 B A d y - U* A N c ) 5{y - y'). (3.22) 

The distinguished role of the Lagrange multiplyers manifests itself in the fact that only 
the a = B component of (|3.11|) forms the first order differential operator while the 
a = i components are ultralocal in y. 

Sometimes a composition of differential operators results in an ultralocal operator. 
Here is one important example that follows from the transformation property of the 
momentum. On the one hand it is given by the canonical transformation 



3 Tangential to the brane components S M | = £ M should not necessarily vanish because these dif- 
fcomorphisms do not shift the boundary. In the canonical formalism this property follows from the 
linearity of the momentum constraints in p, due to which relevant surface terms identically vanish 
for arbitrary S**|. 
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On the other hand it can be obtained by the gauge transformation of the Lagrangian 
expression for the momentum (J3.14)) . The metric variation of the latter has the form 
of the differential operator acting on 5G a (y) = (5g l (y),SN A (y)), 

5G dL»(g,g,N) ^ w s {d)5G a {y) (3 24) 

In essence Wf a (d) here is a part of the Wronskian operator ()2.25j) ) associated with the 
part of the operator (j2.16|) . F ab = F AB ' CD (V) 5(X, X'), without the gauge-breaking 
term. Substituting the gauge variation of the metric one therefore has 

<Pp°(<7, g,N) = Wf a 6 E G a (y) = Wf a R a A E A (3.25) 

and comparing (J3.23|) and (J3.25|) finds that the composition of the two first-order dif- 
ferential operators is ultralocal in y (and contains at most the derivatives with respect 
to brane coordinates) 4 



dg* 



(3.26) 

p=p°{g,g,N) 



The final comment of this section concerns the diffeomorphisms tangential to the 
brane. In the set of R l A = (R 1 ^, R l ± ) their generators R l = R l ^(g) are independent of 
g % and correspond to the momentum constraints linear in p, 

H ll (g,p)=R i tlPl . (3.27) 

By construction they leave the brane part of the action invariant 

Rl — = 0, (3.28) 

which is the brane, d-dimensional, analogue of (|3.3jK 



4. Gauge-fixing procedure 

The relativistic nature of the bulk gauge conditions (|2.1j) which in condensed notations 
we denote by 

H A = H A (G a (y),G a (y)) (4.1) 

4 The equivalence of the canonical and Lagrangian gauge transformations of momenta, 
\d a p)\p°!g,g,N) = (p° (g , g , N)) , holds as it follows from l|3.16|l only up to terms proportional to 
equations of motion ^J]. In what follows the bulk background always satisfies on-shell condition 
which justifies the identification of canonical and Lagrangian versions of gauge transformations. 
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implies that they necessarily depend on the "velocities" of the Lagrange multipliers 
N A = d y N A , so that the following matrix is nondegenerate 

dH A 

4 = -^, deta^O. (4.2) 
dN B 

As a result the corresponding Faddeev-Popov operator in view of the derivative nature 
of the transformations ()3.20|) for N A is of the second order in d y . This follows from the 
definition of this operator 

Q A B = 5 -^R a B= H A a (d) R a B (d) 5(y - y'l (4.3) 

S J^ = H A a =H A a (d)5(y-y'), (4.4) 

where H a is the matrix of linearized gauge conditions (|2.21|1 ). This functional matrix 
in relativistic gauges is the first-order differential operator, H A c (d y ) = —a A d y + 
whence in view of (|3.22|) Q A B (d y ) = —a B dy + ... . Such a second-order operator admits 
zero modes parameterized by their boundary conditions, corresponding to the residual 
transformations discussed in Introduction and in Sect. 2. 

To fix these residual transformations we impose gauge conditions on the brane 
metric (|2.3jh x M (flO = 0- They force these boundary conditions to vanish, provided the 
brane Faddeev-Popov operator, which in condensed canonical notations reads as 

is nondegenerate. As discussed in Introduction, the brane gauge conditions can be 
chosen background-covariant from the c?-dimensional viewpoint, as is the case of the 
brane DeWitt gauge (|1.9|) . However, they are imposed only on g % = g a p{x) and do not 
involve g\ so that they can be considered unitary from the viewpoint of the y-time 
canonical formalism. 

Thus, the overall gauge-fixing factor in the path integral (jl.8|) takes the form 

(gauge-fixing) = 6 [H A (G)} S^g)) M H , X [G], (4.6) 

where we use different condensed notations for the functional delta-functions in the bulk 
(J2.6J) and on the brane (|2.7j) . The measure factor Mh, x [ G ] is determined according 
to the standard Faddeev-Popov procedure by the following functional integral over the 
full diffeomorphism group 

(M X ^ H [G] y 1 = J DZ5[H A (G E )]5(x»(gt)). (4.7) 
3^1=0 
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To calculate it we write the infinitesimally transformed bulk and brane metrics as 
G$ = G a + R a A E A + 0[ E 2 ] and g\ = g i + ^ + 0[ £ 2 ], and then decompose in (|T7jl 
the integration over E A (X) into the integration over the bulk field with fixed boundary 
values E A (X) \ = (0,£ M (x)) and the subsequent integration over 

j DE5[H A {G,)]5{ X "{g^)) = f di5{JX) j DE5[Q A E B ] 

H^l =0 £ x | =0 

3" I _ 

= (detJe)" 1 J DE5[Q A Z B ]=(detJtr 1 (Det D Q A y 1 . (4.8) 



The result is structurally very simple - it factorizes into the product of functional 
determinants of brane and bulk Faddeev-Popov operators, and the latter, DetDQ^, 
is calculated subject to Dirichlet conditions for all components of E A \. This directly 
follows from the gauge-fixing procedure on the brane. 

The Dirichlet type functional determinant is determined by the variational formula 

6 In Det D Q A = Q~ D 1B 6Q A = f dy Q D Y%', y) 5Q A (d y ) , (4.9) 

J y=y 

where SQ B A denotes the variation of the operator under generic change of its coefficients 
and Qd^ i s the Dirichlet Green's function of the ghost operator, defined by 

Q c B (d y )Q D 1B A (y,y') = S%6(y - y'), 

Q D 1B A (y,y')\ y = o. (4.10) 

Thus, this derivation yields one of the main results of this paper (j2.5|) and confirms 
boundary conditions for the ghost propagator in the bulk. 

As a result of transition ()2.8j) - (j2.11j) the quantum effective action takes the form 

e ir = J DG^[G] exp [iSg[G]} det J£ Det D Q^, (4.11) 

S^G] = S B [G] + S^iG] + S h [g] + S* h [g], (4.12) 

with the full gauge-fixed action S'gffG] including the bulk and brane gauge-breaking 
terms which read in canonical condensed notations as 

S%[G] = ~\ J dyH A (G,G)c AB H B (G,G), (4.13) 



s h Ag) = -^{g)c^x v {g)- (4.14) 



1 
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The last comment of this section concerns the local measure denoted in (|4.11j) by 
fi[ G ]. Its contribution is not very important for practical purposes, because its function 
solely consists in the cancellation of strongest power divergences of the path integral. 
Nevertheless, for completeness we briefly discuss it here. 

The local measure is determined by another type of bulk and brane foliation - 
the one with the slices of constant physical time t = X°. This foliation determines 
the chronological ordering in the unitary evolution in physical (not fictitious) time. 
It gives rise to the local measure as a contribution of the gaussian path integral over 
the physical momenta conjugated to temporal velocities 8qG (rather than to d y G). 
As shown in |Appcndix when the brane action has the kinetic term 5 , the full local 
measure factorizes into the product of the bulk and brane measures 

li[G] = iMB[G]ii b (g). (4.15) 

The bulk measure looks as follows. In the foliation of the bulk by spacelike slices 
of constant X° = t, X = (t, x, y) (where x denotes spatial coordinates among brane 
coordinates x), the bulk metric can be decomposed into the corresponding spacelike 
metric Q(t,x,y) and bulk lapse and shift functions N{b,yL,y) ~ GAt(X), Gab(X) ~ 
(Q(t, x, y), J\f(t, x, y)). The bulk measure then reads as 



mb[g]= n 

X,X£h 



, d 2 L B {Q, d t Q,M) , t 
det— —————— detc^ 

d{d t Q) d{d t Q) 



1/2 



(4.16) 



It also absorbs the determinant of the gauge-fixing matrix cab which is generated by the 
transition (|2.8|) - (|2.11j) to nondegenerate gauges and which we also consider ultralocal 
in X. 

Similarly, in the construction of the brane measure one has a foliation of the brane 
by spacelike (d — l)-dimensional surfaces x = (i, x), which leads to ((d — 1) + l)- 
decomposition g(x) = (q(t, x), n(t, x)) into the space metric q(t,~x) and brane lapse 
and shift functions n(t, x). The brane Lagrangian then depends on velocities of only 
the g-metric coefficients, L^(g,dg) = L\,(q,dtq,n), and the brane measure takes the 
form 



M0. 



n 



d 2 L 



det " b detc 
o{d t q) d{d t q) 



1/2 



(4.17) 



5 Which is not always the case - in the Randall-Sundrum model, for example, the brane part of the 
action is kinctically inert, because it contains only the brane tension. 
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5. Bulk wavefunction representation of the brane 
effective action 



As a next step we decompose the full integration in (|4.11j) into the integration over the 
bulk metric Gab(X) subject to fixed induced metric on the brane g a p{x) = G a p(X) \ 
and the subsequent integration over <i-dimensional g a p(x), 

J DG AB {X) = J dg aP (x) J DG AB {X) = J dg J DG. (5.1) 

Gap | =g a p(x) G i \ = g i 

The result looks as the Feynman-DeWitt-Faddeev-Popov functional integral [Tlj 

e lF = J d gf , h [g}ex V i(S h (g) + S* h (g)) detJ»* B (g) (5.2) 

for the purely <i-dimensional system with the brane action S h [g a p] but with the inser- 
tion of the functional \P B (g) which we will call a wavefunction of the bulk spacetime 



& B (g)= J DG/2 B [G] ex P i(S B [G] + S* h [G]) Det D Q^. (5.3) 

G i \=g i 

This function is well known from quantum cosmology as a path-integral representa- 
tion of the solution of the Wheeler-DeWitt equations - quantum Dirac constraints on 
a quantum state in the canonical quantization of gravity. At various levels of rigor and 
in various contexts it was built in a path-integral form in [TBJ 113 EEE1 E3 120] • In partic- 
ular, in the d = 3 context it was obtained in ^H] where the initial and final spacelike 
Cauchy surfaces played the role of branes with specified 3-metrics. The metrics served 
as functional arguments of the two-point kernel of the unitary evolution interpolating 
between these two 3-surfaces. 

In (15 .3j) the role of the boundary is played by timelike brane(s), and the canonical 
formalism in real time is generalized to the case of the "evolution" in the direction 
transversal to the brane. Certainly, no sensible unitarity or causality can be ascribed 
to this evolution. Apparently this construction can be generalized to a purely Euclidean 
case when the full boundary can have various topological and connectedness properties. 

The bulk wavefunction (|5.3J) has the same properties as the cosmological wavefunc- 
tion of PHI 120] ■ First of all, it is independent of the choice of the bulk gauge H A , 



Wb(s) = 0, 



(5.4) 
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though H A explicitly enters its construction. This a typical property of the Faddeev- 
Popov functional integral on shell (or in the absence of sources located in the bulk - 
the only source of \&b{9 ) is the brane metric g = g a p{x) ). A formal proof is based on 
the change of the integration variable in ()5.3|) 

G a -> G a 3 = G a + R\ H A , E A = Q^ A SH B , (5.5) 

simulating the change of gauge conditions H A — > H A + 5H A . Due to Ward identities 
for ghost and gauge Green's functions and gauge invariance of the bulk action, this 
transformation is identical and leads to (|5.4j) . Important point of the transform (J5.5|) is 
that the Green's function of the Faddeev-Popov operator Q~g has Dirichlet boundary 
conditions, and therefore S A | = 0, so that this transform does not shift g l on the brane 
(as it should). It only shifts lapse and shift functions N c \ because of the y-derivative 
in R^(dy) acting on Ii A (y), but this is not dangerous because the Lagrange multiplyers 
N c \ are integrated over at the boundary, rather than being fixed like G l \ = g % . 

Important consequence of this integration over N c \ is that ^(fl 1 ) satisfies the 
analogue of the Wheeler-DeWitt equations 

H A (g,d/idg)*B{g) = 0, (5.6) 

where ff^^, d/idg) are the operators of quantum Dirac (momentum and Hamiltonian) 
constraints. They follow from their classical counterparts by some nontrivial operator 
realization which is formally known only in the linear in H approximation [2*T]. 7 

The momentum components of (|5.7j) . A = fi, describe the gauge invariance of g ) 
under <i-dimensional diffeomorphisms. In condensed notations this reads as 

This property can also be independently proven by the transformation of integration 
variable similar to ()5.5|) but with E^] = C ^ - the diffeomorphism tangential to 

6 This is a point of departure from where the effective action is a functional of N c \ and, therefore, 
depends on the choice of the gauge in the bulk. 

7 It is not obvious, though, that the realization of [2] for the temporal constraints would work for 
brane constraints, because the latter govern an artificial evolution in the y-time which is not related to 
unitarity and causality. In particular, the choice of the local measure in (|5.3|l is not related to y- foliation 
of the bulk and does not yield the scalar density nature of ?Pb ( 9 ) under the reparameterizations of 
g. This choice is tightly related to the operator realization of quantum Dirac constraints. For the 
temporal constraints of |21j the choice of measure guarantees that the cosmological wavefunction is a 
scalar density of 1/2- weight, and the operators Ha (g, d/idg) act covariantly on this density. No such 
properties hold for I^b ( g ) - a scalar in g rather than a density, see |Appcndix A| 
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the brane [THj. This property in its turn guarantees the %-gauge independence of the 
on-shell brane action r, 



5 x r = 0. (5.8) 

Indeed, the <i-dimensional analogue of the change of integration variables (J5.5|) in the 
path integral (|5.2|) gives this property in virtue of gauge invariance of S h ( g ) and \Pb( g ) 
and the Ward identities for the brane Faddeev-Popov Green's function. 



6. Semiclassical expansion 
6.1. Wavefunction of the bulk 

Semiclassical expansion for ^(fiO is based on the stationary point of the gauge- fixed 
action in the bulk. As a result of integration by parts, the variation of the action 
includes the bulk and boundary terms 

6 (S»[G] + SllG] ) = JdX - C DB H B )S(P{X) 



B 



D 



dG a dG a 



(6.1) 



which separately should be equated to zero. Now, take into account that for fixed 
G l \ = g l the boundary variation 8G % \ = and note that Lb is independent of N A and 
the matrices dH D /dN A = -a D A (cf. Eq.flOJ) and cdb are invertible. Therefore, the 
stationarity conditions reduce to 

5S B 5H D uB 

Jghx)' ~Jghx) CdbE =0 ' (6 ' 2) 

H B \ = 0. (6.3) 

Functionally contracting ()6.2j) with R a A and using gauge invariance of the bulk action 
(J3.3|) we get for the gauge conditions the homogeneous equation with the second-order 
Faddeev-Popov operator 

Q^(d)c DB H B (X) = 0. (6.4) 

Since in view of ()6.3|) they have Dirichlet boundary conditions, the solution is identically 
zero everywhere in the bulk. Thus, eventually the stationary point of the bulk action 

G a (g) = G AB (X)[g a p(x)} (6.5) 
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is a solution of classical equations of motion in the bulk gauge with fixed brane metric 
at the boundary, as advocated in Eqs. (J2.13)) - (J2.15j) of Sect. 2. 

The quadratic part of the action on this background - the variation of (|6.1|) - reads 

as 

^5 2 (S B [G] + S%[g]) = l -5G a (F ab 5G b ) + i 5G a W ab (d) 5G b \ , (6.6) 
where F a b is the condensed notation for the operator (j2.16|) 

F a b = S a b — Ha c AbH^, (6.7) 



5 2 S B [G] 
5G*(y)5G b (y>) 



Sab = Sab (d y ) 5{y - y') = s „ n/ ^ (61 



and Wab (9) is the corresponding Wronskian operator (|2.25jl defined by the variational 
relation 

W ab (d) 5G b (y) =s(?i-- ^C- cabH 1 
1 ; KU> \dG a dG a 

= (W! b (d) -^2 CAB H B b (d)) SG\y). (6.9) 



H A =0 

dH A 



Together with the Wronskian operator of the invariant action W^ b (d), see Eq. (13.24)1 . 
it includes the contribution of the gauge-breaking term. Wab (d) participates in the 
Wronskian relation (|2.25)1 and also satisfies an additional relation that can be obtained 
by a single integration by parts of the derivatives in F a b (9) 

ttFabrt = 0lFa602 ~ <f>l Wab (0) $ | (6.10) 

Here F a b implies that the derivatives in F a b (d) are integrated by parts one time to 
form the bilinear combinations of the first-order derivatives of 0^2 (for 2 = (f) 1 this is 
just the Lagrangian quadratic in <fi and dip). With this relation the quadratic form of 
the action (|6.6j) takes the form 

1 -5 2 {S B [G] + Sl[g]) = l -5G a F ab 5G b . (6.11) 

The one-loop contribution of this form to the path integral ()5.3j) is the gaussian 
functional integral over the perturbations 5G a = 5Gab(X) subject to Dirichlet bound- 
ary conditions for the brane metric components 5G l \ = 0, while 5N A \ are integrated 
over in the infinite limits. As shown in Appendix B this integral yields the functional 
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determinant of F a t> subject to mixed Dirichlet-Neumann boundary conditions, and the 
one-loop wavefunction of the bulk has the final form 



'Mtf) = -J72 exp US [G 

(Det DN F ab ) 1/2 V 



(6.12) 

G = G (g) 

Here the Dirichlet-Neumann functional determinant DetoN-^ab is determined by the 
variational relation 

6 In Det DN F ab = 6 Tr DN ln F ab = 6F ab G b ^. (6.13) 

The notation 5F ab means arbitrary variations of the coefficients of the operator, and, 
similarly to Eg. (jfi.lOj) . the double arrow implies symmetric action of two first-order 
derivatives of SF a b = 5F AB,CD (V) on both arguments of the Green's function G*d N = 
AB (X,X') (before taking the coincidence limit X' = X implicit in the functional 
trace operation). 

The most peculiar element of the variational definition (J6.13j) is G b ^ N - the Green's 
function of the operator F a b(d) subject to the mixed set of Dirichlet and generalized 
Neumann boundary conditions (|2.18J) - ()2.20|) . In canonical condensed notations this 
boundary value problem reads as 

F ca (d)G%(y,y')=5 b c 5(y-y'), (6.14) 
G% N (y,y')\ y = 0, (6.15) 

H A a (d)G a Uy,y')\ y = o. (6.16) 

Thus, contrary to the purely Dirichlet ghost propagator the metric propagator has 
Dirichlet conditions only for i = (a(3, X) components. The rest of boundary conditions 
()6.16|) belong to the generalized Neumann type. This is a consequence of the fact that 
lapse and shift functions in the path integral for ( 9 ) are integrated out on the brane. 

6.2. Brane effective action 

Now substitute one-loop bulk wavefunction (|6.12|) into the path integral (|5.2j) and 

calculate it by the stationary-phase method. The stationary point of the overall tree- 
level phase satisfies the following equation 



|-(5 B [G (,)] + 5 b (,) + 5 g b b (^)) = |^, 



dS h <9y m 



because S B [Go(g) ] plays the role of the Hamilton- Jacobi function in the y-time canon- 
ical formalism, and its gradient in g % yields the canonical momentum on the brane 
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p1(G, G) = 8L B /dG l 



A 5 B [G ( )]= ^B 



(6.18) 



In essence ()6.17|) represents the generalized Israel matching condition relating the 
extrinsic curvature of the brane, K af3 (x) ~ dL B /dG a/ 3(x), to the brane stress tensor 
dSb/dg 1 = SSb/Sg a /3(x) in the presence of the brane gauge-breaking term. To handle 
the latter, contract (|6.17|) with R l and take into account that in view of (J3.18|) and 
(|3.27j) the first term in the contraction vanishes on bulk shell, 



R] 



, dL, 



5S l 



dG l 



5N» 



= 0, (6.19) 

G=G (g) 

whereas the second term vanishes in view of rf-dimensional covariance of the brane 
action (|3.28j) . so that finally 



p« _ ^ _ -,0 

M Qgi ^ X 



J x a c XaX a = 0, 



(6.20) 



whence x M = due to invertibility of J x and c\ a . Therefore, we get the following 
system of equations for the stationary point go in the <i-dimensional brane gauge 



dL 



B 



dG l 



+ 



dS l 
dg i 



0. 



9 = 90 



0. 



(6.21) 
(6.22) 



To find the quadratic form of the tree-level action on the background of go we 
would need the derivative of G^{g) with respect to g l . This quantity satisfies the 
linearized version of the boundary- value problem (j2.13|) - (j2.15|) . In condensed notations 
this problem for the bulk perturbation <f) a ((p) induced by the perturbation if 1 on the 
boundary reads as 



Sa b (d)<j> b (y) 
Hi(d) <f> b (y) 



0, 
0. 



(6.23) 
(6.24) 



To solve it notice that the Ab component of the Wronskian operator (|6.14j) is entirely 
determined by the gauge-breaking term of the full action 



W Ab (d) 



dH 



B 



dN A 



c BD H u b {d)=a AD H u b (d) 



a A °BD- 



(6.25) 
(6.26) 
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Therefore the problem ()6.23j) - ()f).24)) can be identically rewritten in the form 

F ab (d)4> b (y) = 0, (6.27) 

W Ab (d)<P b \=0, <^|=^\ (6.28) 

where the nondegenerate operator F a b replaces the degenerate operator Sab defined 
by (|6.8j) . This is a problem with the (inhomogeneous) Dirichlet-Neumann boundary 
conditions. By using the Wronskian relation (|2.26|) its solution can be presented 
in terms of the Dirichlet-Neumann Green's function of the above type 

r D M = -G%(y, y') w u (£>') | y = -g% w u \ v\ (6.29) 

whence 



Using (J6.3U)) together with (|6.18|) and (|3.24|) one immediately obtains 

^[f^ =-W*aG%W bk \\, (6.31) 

where the double vertical bar denotes the restriction to the brane with respect to 
the both points of the two-point kernel, cf. Eq. (J2.24j) of Sect. 2. This restriction, 
in particular, allows one to replace W f a by Wia in view of (j6.16J) . As a result, the 
quadratic part of the full tree- level action on the background go, g = go + Sg, takes the 
form 

\ 51 {S B [G (g) ] + S b (g)+S h gh (g)) = 1 -6g* 5g k , (6.32) 

where F™ is the full brane-to-brane operator introduced in (|2.23|) 

FT = - W ia G^n W bk 1 1 + 5& - X? c^xl (6.33) 

with the gauge-fixed contribution of the brane action ()2.27|) 

b b d 2 S h dx 41 

K ik = S ik - Xi c^Xk, s ik = Q g iQ g k ' & = ~QgJ- ( 6 - 34 ) 

Finally, substituting ()6.12|) into ()5.2|) and taking the Gaussian integral over Sg l in 
the vicinity of the stationary point go one finds in the one-loop approximation 

ir(q ) /fg Det DQB f*b det J£ ■/ C Br^i , oW \\ 

6 = f Det F V' 2 (deiF™V /2 eW%{ [ ] ig0)) ' ( ] 

(Uet DN ^ a6 j [det* ik ) G=G ( go ) 

The preexponential factor here confirms the decomposition property for the one-loop 
effective action ()2.22|) advocated in Sect. 2 (modulo 5(0)-type terms generated by the 
local measure (|4.15j) ). 
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7. Ward identities and gauge independence 



It is important that the both bulk and brane preexponential factors in (|6.35j) are 
separately gauge independent. For the bulk prefactor this is a direct corollary of the 
gauge independence of \P-Q(g), (J5.4j) . However, it is worth showing explicitly how this 
property works in virtue of the Ward identities for the bulk propagators. 

Bulk Ward identities (for tree- level propagators) arise as follows [20] • Functional 
differentiation of (J3.3J) shows that on shell, that is on the background satisfying clas- 
sical equations of motion, the functional matrix S ab is degenerate because it has zero- 
eigenvalue eigenvectors - the gauge generators 

5R a 

As a consequence the operator F ab satisfies the relation 

RjiFab = —Qa B c BcH c b (7.2) 

which can be functionally contracted with matrices of the gauge G bc and ghost Q^ 1 A 
Green's functions. Integration by parts of the derivatives in the operators R a A = R a A (d) 
and Q A = Qa(9) does not produce additional surface terms in view of the boundary 
conditions for the propagators. Thus, one arrives at the identity relating the Dirichlet- 
Neumann Green's function of the metric operator to the Dirichlet Green's function of 
the Faddeev-Popov operator 

c AB H B b (d)G b Uy,y') = -Q- A 1B {y,y')R£{d y ,). (7.3) 

It is important that the Dirichlet conditions at y belonging to the boundary on the 
right hand side match with the Neumann conditions (|6.16j) on the left hand side. 
Therefore, using the variational definitions of the functional determinants one has 

6 H In Det DN F ab = 5 H F ab G% = -2 c AB H B G% 5H£, (7.4) 
S H \nDet D Qi = Q A lB 5 H Qi = Q A 1B R B b SH b A , (7.5) 

so that in virtue of ()7.3|) the bulk part of the one-loop effective action is if-independent 



|Tr DN lnF + Tr D InQ 



Bulk gauge obviously participates in the construction of the brane part of the effec- 
tive action, but the latter turns out to be also independent of H A . To show this, note 
first of all that in view of (IfHfij) W ia G a ^ N Wbk\ \ = W S ia G a £ N W S bk \ \- This replacement 
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of the Wronskian operators W(d) by their if -independent analogues IV s (<9), cf. Eq. 
(HOD , implies that 



6 H F™ = - Wi (5 H G%) W s bk 1 1 . (7.7) 
The variation of the Green's function accounting for its boundary conditions is derived 



in 



Appendix C Using (jC.5|) one has in view of ()7.3|) 



r t^iab ^yac r j-i r~idb _ c\r~i{o-c jj A r ttB /^db) 

Offl-f DN ~~ — ^ DN °H-r cd <-r DN — /Or DN -fJ c CAB On d U DN 

= -2S5Q-^c5>,G* ) N , {7.1 
so that in virtue of the identity ()3.26j) the argument of the ghost propagator in 



8 H F™ = 2wLR a AQ- B A 5H B d G%W bk) 



(7.9) 



is not differentiated at the boundary because dHA/dq 1 is ultralocal in y, and the 
whole variation vanishes because of the Dirichlet boundary conditions for the ghost 
propag ator, (dHA/dq*) Q'^ A \ . = 0. Finally, 

5 H (~ tr In F DN + tr In J^j = 0. (7.10) 

Let us turn to the dependence of the effective action on the brane gauge As 
in usual gauge theories, the effective action is gauge-independent on shell, that is with 
the sources switched on the brane off. 8 To see this in the one-loop order, note that the 
tree-level brane action is identically invariant with respect to <i-dimensional (brane) 
diffeomorphisms 

^^-( sB l G o(9)] + S\g))=0. (7.11) 
The differentiation of this identity at g = go (mass shell (j6.21|0 then gives 
ijj, (- W ia G% W bk 1 1 + Si ) = 4 (S B [ G (g) ] + S\g) 



9 = 90 



d 

Qgk Qgl 



dgidg 1 

(s B [G (g)] + S b (g)) =0. (7.12) 

'9 = 90 



8 By construction, there are no sources in the bulk, so that both S'bC?) and P are universally 
independent of the bulk gauge H, because the latter is not sensitive to the inclusion of boundary/brane 
sources. On the contrary, the brane gauge-fixing procedure is vulnerable for gauge non-invariant 
sources on the brane. 



27 



This means that the the brane-to-brane operator (jfi.HHJ) is non-degenerate entirely due 
to its brane gauge-breaking term. Therefore, the following relation holds between the 
operators and J 1 ^ 

R^Ffk 1 = -J V n c va Xk- ( 7 - 13 ) 

Consequently, their Green's functions satisfy the brane Ward identity 

J'l* R\ = -c vaX tG k \ (7.14) 
F? k N G km = ap, (7.15) 

where G kl is the brane-to-brane propagator. The relation ()7.14|) is a direct analogue of 
the bulk Ward identity ()7.3|) . except that we don't have to care about integrations by 
parts in x, because the relevant surface terms are vanishing at the brane infinity x° —> 
±00 (or missing in Euclidean context in view of the closed nature of the boundary). 
Thus, using the variations 

6 x F? k * = -25x1 c - Xt), S X J; = H 6tf (7.16) 

we immediately obtain in view of ()7.14j) the on-shell gauge independence of the brane 
effective action 

5 X M tr In F DN + tr In jj = b X \ c» a xt G kt + 5£ J" 1 " = 0. (7.17) 

8. Reduction of the Dirichlet-Neumann problem to 
the Dirichlet type 

If we would apply the stationary phase method directly to the path integral (j4.11j) 
without using the bulk wavefunction representation (|5.2|) . then the one-loop contribu- 
tion of the metric field would be determined by the functional determinant subject 
to Neumann boundary conditions for all metric components. So the decomposition of 
the integration procedure (|5.1|) exercises a half-way reduction to the Dirichlet problem, 
because the resulting boundary conditions are of a mixed type. As it was discussed 
in Introduction and in the construction of the Dirichlet Green's function is much 
easier than the Neumann one, so actually no calculational advantages are gained when 
this reduction is incomplete. Therefore, it is worth making a further reduction from 
the Dirichlet-Neumann problem to the purely Dirichlet one. 

For this purpose introduce the Gaussian path integral over metric perturbations 
5G a = (j) a = (f> A ) with slightly more general than in (J6.11|) - inhomogeneous - 
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boundary conditions, = ip l 



Z(<p*) = j £>0exp(-S[0]), (8.1) 

<j>'\ = p' 

S[</>] = ± J dy^F ah <$>. (8.2) 

Its stationary point ^dnIv 9 *) is determined by the Dirichlet-Neumann problem with the 
inhomogeneous boundary conditions (|6.27|) - (j6.28j) and has the form (J6.29|) . As shown 
Appendix B[ the gaussian integral itself equals 



in 



Ztf) = (DetG DN ) 1/2 exp(-S[0 DN (^)]) ! (8.3) 
where the action at the stationary point is 

S[<fcxtf)\ = - l -ip l (WiaG&WnW) </\ (8.4) 

Alternatively it can be calculated by integrating first over the fields subject to 
Dirichlet boundary conditions for all <p a with the integration over ip A = (p A \ reserved 
to the last 

Ztf) = 1 J] d f A J D<P exp(-S[0]). (8.5) 

A <f> a \=tp°- 

The inner integral 

J D^expi-S^}) = (DetG D ) 1/2 exp(-5[0 D (^)]) (8.6) 

is determined by the Dirichlet problem for the field <$h(<p), -Fa&^D = 0, 0^ | = (p a , 
having as a solution the following expression in terms of the Dirichlet Green's function 
Gd of the operator F ab 

r D ^) = -Ggw bc <p c \ c . (8.7) 
The action in ()8.6|) equals 

S[M^)] = \^ a D ab ^\ (8.8) 

D ab = -W ac GgW db \\, (8.9) 

so that the gaussian integral over the boundary field ip A in (|8.5j) is saturated by the 
saddle point ip A of the above quadratic form in ip A , (p A = —G AB DBi<p\ where Db% is 
the .Bi-block of the operator ()8.9|) and G AB is the inverse of its v4I?-block 

D AB G BC = 5 C A (8.10) 
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(note that the operator (J8.9)) differs from the kernel of the action ()8.4j) by the type of 
the Green's function - Dirichlet vs Dirichlet-Neumann one). 
Thus, integration over tp A in (J8.5J) gives 

Z{^) = (DetG D ) 1/2 (det^)~ 1/2 exp(-5[0 D (v9 O )]), (8.11) 

S[M<Po) \ = \<P i - D ^ GAB D Bk) <p k - (8.12) 
Comparison with (|8.3j) - (|8.4|) then gives the following two relations 

-( WG DN W\\). k = D ik - D iA G AB D Bk , (8.13) 
/ \ 1/2 / \ 1/2 / \~ 1/2 

( Det Gdn J = ( DetG D J (detD AB ) . (8.14) 

They allow one to completely reduce the algorithm ()6.35|) to that of the Dirichlet 
boundary conditions. In particular, the substitution of ()8.13|) into the expression for 
the brane-to-brane operator , (|fi.33|) . expresses its functional determinant as 

det F° N = det ({D lk + n tk ) - D iA G AB D Bk ) = det F° ( det D AB ^j ~\ (8.15) 

where F® b is the brane-to-brane operator of another type - it acts on the space of all 
metric perturbations on the brane 5G a \ = (5g\SN A )\, a = (i,A), and is built of the 
Dirichlet Green's function rather than the mixed-type one 

= -(WG D W\\) ab + K lk 5l5 k b . (8.16) 

Its ifc-block is given by Eq. ()2.29|) in Sect. 2. 

In view of the relation (|8.15|) the one-loop contribution of the bulk and brane 
gravitons in (|6.35|) can be decomposed into the product of other two bulk and brane 
factors which are entirely based on the Dirichlet-type objects 

/ \ -1/2 / \-l/2 / \ -1/2 / \ -1/2 

(Det DN Fj (detF° N J =(Det D Fj (detF°J . (8.17) 

In fact, this decomposition literally repeats the Neumann-Dirichlet reduction suggested 
in p. The corresponding decomposition of the one-loop effective action has the form 
(J2.28J) presented in Sect. 2. From calculation viewpoint it is simpler than ()2.22j) . but it 
destroys manifest gauge independence of the combined bulk-brane diagrammatic tech- 
nique. Whereas in (|2.22j) both bulk and brane parts are separately gauge- independent, 
no such property holds for ()2.28J) - only the sum of bulk and brane terms is independent 
of the bulk gauge conditions. 



ab 



Di k + ^ik D iB 
D A k D AB 
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9. Conclusions 



Thus we have an exhaustive formulation of the Feynman-DeWitt-Faddeev-Popov gauge- 
fixing procedure in gravitational systems with branes or boundaries. This procedure 
incorporates a special choice of gauge conditions in the bulk and on the brane, which 
separately fix the bulk and brane diffeomorphisms. Also it establishes the boundary 
conditions for the corresponding ghost factors and allows one to construct a special 
bulk wavefunction representation of the brane effective action (|5.2j) - ([5.3|) . The bulk 
wavefunction satisfies the generalized Wheeler-DeWitt equations with respect to the 
induced metric of the brane (|5.7|) . They might, perhaps, serve as a basis for non- 
perturbative methods in brane models, alternative to semiclassical expansion [T%] . 

We derived the boundary conditions for propagators of the Feynman diagrammatic 
technique in brane models and explicitly built the brane effective action in the one- 
loop approximation. Similarly to non-gauge (quantized matter) models with branes, 
considered in |]Q , this one- loop action can be decomposed into the sum of bulk and brane 
contributions both in the graviton and the ghost sectors (|2.28j) . The bulk contribution 
is represented by the usual functional determinant of the propagator in the (d + 1)- 
dimensional bulk, subject to Dirichlet boundary conditions on the brane. The brane 
contribution is given by a similar determinant of the brane-to-brane operator ()2.29|) in 
the surface c?-dimensional theory. In the graviton sector this operator has a nonlocal 
(pseudodifferential) nature, whereas in the ghost sector for local gauges this is always a 
local Faddeev-Popov operator associated with gauging away the d- dimensional surface 
diffeomorphisms. This property follows from the generic bulk-brane factorization of 
a gauge-fixing factor (J2.5j) in the path integral for gravitational brane models, which 
holds beyond loop expansion. 

Linear algebra manipulations and Gaussian integrations which underly the above 
results look innocent at the calculational level adopted in this paper. The obtained 
algorithms are, however, marred by ultraviolet divergences and should be regulated 
by some covariant technique. The efficiency and correctness of these calculations was, 
nevertheless, confirmed within the dimensional regularization in pQ. In particular, 
correct expressions for surface Schwinger-DeWitt coefficients were obtained in pQ and, 
thus, guaranteed correct renormalization of logarithmic divergences in simplest brane 
models. 

Such calculations should also apply in gravitational systems considered above, and 
this is a subject of forthcoming papers. We plan to make a synthesis of the Neumann- 
Dirichlet reduction with a systematic curvature expansion method in brane models in 
order to alleviate their formalism to the level of universality of the Schwinger-DeWitt 
technique |2*3*1 EM] . This is important for various applications in quantum gravity and 
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cosmology, including quantum consistency of brane models [H] , the conformal anomaly 
mechanism of the dark energy, that might be facilitated within the brane concept of 
extra dimensions [7j, and the others. 



Appendix A The local measure in brane models 



Let the bulk and brane terms of the action both have kinetic terms quadratic in ve- 
locities with coefficients a(X) and b(x). These terms can be rewritten as one integral 
over the bulk 



(A.l) 



dXa(X) (d t (f)(X)) 2 + I dxb(x)(d t ip(x)) 2 

dX dX'G(X, X') d t (f)(X) d t (f)(X') 

with the overall coefficient of the quadratic form in field velocities 
g(X, X') = 8(x, x') [a(X) 8(y - y') + b(x) 8{y) 8{y') ] . 



(A.2) 



Here <j)(X) and <p(x) = <f>(x,0) = 4>{X)\ denote the bulk field and its boundary value 
at the brane located at y — 0. 

The full local canonical measure reads in terms of the functional determinant of 
this matrix as 



fj,[a(X), b(x) 



DetG(X, X'] 



1/2 



jjL B [a(X)}iJL b [b(x)] 



(A.3) 



and as we show below factorizes into the product of the corresponding bulk fi B [a(X) 
and brane fjbb[b(x)] measures given by 



Det {d+1) a(X)5(X,X')\ 1/2 
Det(d) a(x, 0) S(x, x') 



and 



fi b [b(x)] = (Det(d) b(x) 5{x,x') ) 



1/2 



(A.4) 



(A.5) 



This factorization follows from the expression for the inverse of the functional matrix 



(E3 



Q-\X % X') = 8(x,x') 



a(X) 



S(y - y') 



b{x) 8{y) 8{y>) 
a(x,0) +8(0)b(x) 



(A.6) 
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When substituted into the variational equation for the logarithm of the functional 
determinant in (jA.3|) this expression yields the decomposition into the following sum 
of two terms 



8 In Det Q{X, X') = J dX dX' 8Q(X, X') Q~ X (X\ X) 
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8 {d+l) (0) J dX In a(X) + 8 {d) (0) J dx In (\ + 5(0)- 



b{x) 



x,0) 



(A.7) 



With the obvious limit 9 8(0) — > oo this sum gives the factorization of the functional 
determinant 



Det G(X, X') = ,\ nVr^Ln Det W b ( x ) 5 ( x > x ')> ( A - 8 ) 



Det (d+1) a(X)8(X,X' ) 
Det(d) a(x, 0) 8(x, x') 

which implies the factorization l|A.3j) - (|A.5j) . In its turn this confirms the factorization of 
the local measure in the brane model (|4.15jl - ([4.17jl . The division by Det^ a(x, 0) 8(x, x') 
in the bulk measure (|A.4|) corresponds to the absence of local factors belonging to the 
brane in the product of (J4.16|) . In the notations of this Appendix this means the 
following chain of relations 

II °w = II a ^y^ = U^n-)H a ^y^ 

= TT — U- exp ( 8(0) [ dy In a(x, y)) = ^^^^f ^ . (A.9) 
lla(x,0) P V J y >o J Bet (d) a(x,0)8(x,x') K } 

Note that this division is responsible for zero density weight of the bulk wavefunction 
\P B (g), defined by the path integral (|5.3|) (see footnote in Sect. 5). 



Appendix B The Gaussian functional integral with 

mixed boundary conditions 

Feynman's calculation [22] of the gaussian functional integral with the action (j6.11|) 
for the case when only a part of integration fields is fixed at the boundary, 

SG a = 0° = (cj)\ cp% 0*1=0, (B.l) 



9 Any regularization implies this limit except the dimensional rcgularization in which 5(0) = 0, but 
in the dimensional regularization the contribution of the local measure is identically zero. 
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can be based on the integral 

Z[F,J]= J D0exp(-S[0,J]), (B.2) 

4 | =o 

S[<I>,J] = f dy ^ 1 -r(y)F ab (d)<p b (y) + J a (y)<f ) a (y)y (B.3) 

where the action is modified by the source term in the bulk. To find the dependence 
of (jB.2|) on J a consider the stationary point 0dn = </>dn[^] of this action with respect 
to variations of <p a satisfying (|B.1J) . 

5S[<t>,J\ = J dy5<p a (F<P + J) a + 5<P A (W<P) A \=Q (B.4) 

Thus, 0dn satisfies the problem with mixed Dirichlet-Neumann boundary conditions 

F ab (f> b DN + J a = 0, (B.5) 

0dn| = 0, H>Aa0D N |=O, (B.6) 
and therefore can be represented in terms of the corresponding Green's function of 

Eqs. dnmD-dnnsi) 

= G&Jb = J dy'G^y') J b (y') (B.7) 



Now make the shift of the integration variable in (jB.2|) by 0dn 5 4> = </>dn + A. 
Under this replacement the action decomposes in the part S[ A, 0] quadratic in A and 
the part independent of A. Linear in A term is absent (both in the bulk and on the 
boundary) in view of the stationarity of the action at 0dn? so that 

S[hJ] = S[A,0] + S[<fo N ,J], (B.8) 

S[(hN,J] = \j a G^J b . (B.9) 

Therefore 

Z[F,J] = Z[F,0] exp(-S[<f> m ,J]). (B.10) 

To find the prefactor, consider the variation of the integral (|B.2|) at J = with 
respect to the operator F ab and make the following set of obvious identical transforma- 
tions using the above equations ()B.9|) and (jB.lOJl 

6 F Z[F,0] = - f D<\> (^^SFabAexpi-Sl&O]) 



hikw/^u) zif ' j] 



J=0 



- l - 5F ab G b ^Z[F,0}. (B.ll) 
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Here 5F a b means arbitrary variations of the coefficients of the operator, and the dou- 
ble arrow implies symmetric action of two first-order derivatives of F a b{d) on both 
arguments of similar to Eq. (|B.3|) . 
Therefore, one gets 

6 F hx Z[F,0] = - X - 5F ab G<g N = -i<51nDet DN F = ^lnDetG DN . (B.12) 

This justifies the variational definition (JbM3|) of the Dirichlet-Neumann functional de- 
terminant of F a b originating from the Gaussian integration with mixed boundary con- 
ditions. 



Appendix C The variation of G^ N 

In view of (|6.25|) H f = a~ 1AC Web, and the boundary conditions in the problem 
(jfi.l4j) - (jfi.lfij) for Gdn can be rewritten in terms of the Wronskian operator Web- Then 
the variation of this problem gives the set of inhomogeneous equations for SG^ 

Fab ^Cp N = — 5 F a b C^ N , 

W Ab <5^DN I a = ~ SWAb Gdn\ A 7 

SG- N \ i = (C.l) 

This set of equations forms the boundary value problem, analogous to ([6.27)1 - ()6.29[) . 
but with the inhomogeneous subset of Neumann boundary conditions 

F(j) = J, 

(W ( f ) ) A \=w A , 0*|=O. (C.2) 



It has a solution 



dX'G a ^(X,X')Jb(X')+ [ dxG a g N (X,e(x))w A (x) 

Jh 

= G"dn Jb + wa (C.3) 

A 

which, when applied to ([C.l)) . gives 

ex-fab fyac c 77 ^ydb syaA XTT/ s~idb ( r~\ \\ 

"^DN _ — DN v"cd t^DN — ^ DN 6 VV Ad <-^DN l^- 4 J 

A 

Therefore, in view of the relation (jfi.K)j) . the variation of the Dirichlet-Neumann Green's 
function is given by the equation 

^dn = ~ G a vj<i$F cd GnN, (C.5) 
used in Eq.(JHJ). 
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